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Abstract: In the MSSM with complex parameters loop corrections to the decays H + — > tb 
and H~ — > tb with t — > bW and W — > Iv lead to CP- violating asymmetries: a decay rate 
asymmetry, a forward-backward asymmetry and an energy asymmetry. We derive explicit 
formulas for them and perform a detailed numerical analysis. We study the dependence 
on the parameters and the phases involved. In particular, the influence of the running 
Yukawa coupling is taken into account. The decay rate asymmetry can go up to 25%, the 
forward-backward and the energy asymmetry up to 10%. 
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1. Introduction 

It is well known that supersymmetric models contain new sources of CP violation if the 
parameters are complex. In the Minimal Supersymmetric Standard Model (MSSM) the 
U(l) and SU(2) gaugino mass parameters M\ and M2, the higgsino mass parameter /x, as 
well as the trilinear couplings Af (corresponding to a fermion /) may be complex. (Usually, 
M2 is made real by redefining the fields.) Non- vanishing phases of these parameters cause 
CP-violating effects. While the phase of ji may be small for a supersymmetric particle 
spectrum of 0(100 GeV) due to the experimental upper bounds of the electric dipole 
moments (EDMs) of electron and neutron, the trilinear couplings of the third generation 
At,b,T are n °t so much constrained and can lead to significant CP-violation [1, 2]. 

In the following, we study CP violation in the decays of the charged Higgs bosons H ± 
within the MSSM. There are three possible decays of H ± into ordinary particles: H + into 
tb , tv and WhP and the CP conjugated ones, where h° is the lightest neutral Higgs boson 
of the MSSM. At tree level the partial decay widths of H + and H~ are equal because of CP 
invariance of the Higgs potential. However, including loop corrections with intermediate 
SUSY-particles, they become different due to the CP violation induced by the complex 
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phases of the MSSM parameters, essentially of A t & r . Quite generally, these phases affect 
the whole Higgs sector of the MSSM substantially [3, 4]. 

A full one-loop calculation within the MSSM was done of the decays mentioned [5, 6, 
7, 8], and the CP-violating decay rate asymmetry 5 = [T H + — T H -]/[T H + + T H -] for 
these decays was calculated. In the case of H + — > tb and H~ —* tb this asymmetry can go 
up to ~ 25%. 

In this paper, we go a step further by including the decay product particles of the top 
quark, see Fig. 1, 

H + ->bt -^bb'W + , 



bt -> bb'W 



XI) 



and 



H + -►&*-> bb'W + -► bb'l + v h , 

H~ ->bi->bb'W- -tbbTPi. (1.2) 

We only consider CP violation induced by the loop diagrams of H^tb vertex. We neglect 
CP violation in the tWb' vertex. In the Standard model it is very small and in the MSSM 
for m H + > mt all two-body decays of top into SUSY partners are excluded kinematically. 





b b 
Figure 1: The Feynman graphs of the processes we study. 

In particular, we will exploit the polarization of the top quark. The top-quark de- 
cays before forming a bound state due to its large mass, so that the polarization can be 
measured by the angular distributions of its decay products. The polarization is very 
sensitive to CP violation. We will consider suitable CP violating forward-backward and 
energy asymmetries by using angular or energy distributions of the decay particles. The 
asymmetries also depend on the sensitivity of b' in (1.1) or the final lepton in (1.2) to 
the top-quark polarization. Further we make a numerical analysis for different values of 
the MSSM parameters. 

This paper is organized in the following order. In Section 2 we present the formalism 
we use. There are subsections devoted to the polarization of the top quark and the CP- 
violating asymmetries. Sections 3 and 4 contain the angular and energy distributions and 
the analytic results for the asymmetries of (1.1) and (1.2), respectively. The numerical 
results are discussed in Section 5. Section 6 contains the conclusions. In Appendix A the 
formulas used for running Yukawa couplings and ht are given. In Appendix B we point 
out an error made in an equation of [5]. 
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2. Formalism 



In order to obtain the analytic expressions for the differential partial decay rates of (1.1) 
and (1.2), we follow the formalism of [9]. In accordance with it, for both of the processes 
we write 

dr± = dr H± drf jt --^-, f = b',i. (2.1) 

E t 1 is the energy of i(i)-quark, and Tt is the total decay width of the t-quark. dTu± is the 
differential partial decay rate of the process — > tb when CP-violation is included: 

dT H± = ^-\M H± \ 2 d<S> H± . (2.2) 

where d& H ± is the relevant phase space element and 

M H+ = iu( Pt )[Y b + P R + Y+PlH-pb) , (2.3) 

M H - = iu(p h )[Y-p R + Y-p L ]u{- Pl ) , (2.4) 
Y t ± =y t + SY t ± , Y h ± = y b + 5Y b ± . (2.5) 

Here y t and y b are the DR running couplings, see Appendix A, 5Y^ b are the SUSY-induced 
loop corrections, which most generally have CP-invariant and CP-violating parts: 

5Y± = SY^ ± \SY% P ■ (2.6) 

dT{t is the differential partial rate of the process t — ► ^/VF 1 * 1 or i ^ 6 / Z ± ^ when the 
t-quark is polarized and its polarization is determined in the former process — > tfr: 



n , (&,tP/) 

1 ± ctfmt-, r- 

(Ptp/) 



is the polarization vector of the t-quark and d&. + are the phase space elements. The 
index / stands for the corresponding fermion (/ = b',l) and aj determines its sensitivity 
to the polarization of the t-quark: 



mf — 2mm 

"' = fef- a ' = 1 - {2i 



In the kinematics of both of the processes (1.1) and (1.2) we work in the approximation 

m 1/ m w — m b/ m t — m b/ m W — > ( 2 - 9 ) 
but we keep m^ ^ in the couplings, where it is multiplied by tan (3. 
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2.1 The t-quark polarization vector 

The polarization four-vectors £f and £° for the considered processes are covariantly given 
by the expressions [9]: 

r = ( a ap _ Pf£\ Tr[M H+ (-A(-p 5 ))M £ +A( Pt ) 7/375 ] 
* V m? / Tr[M i/+ (-l)A(-^)M if+ A(p t )] ' 1 ' j 

r « _ f g «fl PfP? ^ ?M M h - (-A(-Pt))Mg- A(p & ) 7/37 5] f2U , 
' V < ; Tr[M H -(-l)A(-p f )M H -A(p b )] ' 1 ' j 

where 

M H + = Y+P R + Y+P L , M H - = YfP R + r-F L , (2.12) 
M = 7o Mt 7o , A(pt) = j*t + mt . (2.13) 

Thus we obtain: 

$ f = m^Qfc, Qt- b = p% b - te^jfy . (2.14) 

"H 

Notice that the four-vector Q bb is the only four-vector in — > tb that can be constructed 
so that it satisfies the orthogonal condition (£tPt) = 0. The polarization vectors (2.14) have 
CP-invariant and CP-violating parts, and the CP-violating parts are only due to the loop 
corrections: 

<p± = ±<pinv + V CP ^ ^ 2>15 j 

■pinv _ Vt ~ Vb \Q\ 

(Vt + vDiPtPl) - 2m t m b y t y b ' 

v cp = y t Re(SY t CP )-y b M8Y b CP ) ^ 

(% 2 + VDiPtPl) - 2m t m b y t y b ' 
The explicit forms of the individual contributions to Re(5Y t c b p ) are taken from [5]. 

2.2 CP-violating asymmetries 

The CP-violating decay rate asymmetry 5 CP is given by the expression 

6 cp = r+-r_ 

r+ + r_ v ; 

In (2.18) r± are the partial decay widths of H^. 

Next, we construct a CP-violating forward-backward (FB) asymmetry AA CP from the FB 
asymmetries using the angular distributions of the processes 



AA CP = A PB - A PB , (2.19) 



where 



^ B = ifTrf' (2 ' 2()) 



r£= / -dcos9 and r| = / -dcosfl, (2.21) 

J dcosO Jk a cos 6 
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i.e. are the number of particles /antiparticles measured in the forward direction of the 
decaying t/t quarks, etc. 

Analogously, a CP-violating energy asymmetry AR CP can be defined, using the energy 
distributions of the processes 

AR CP = R + -R^, (2.22) 

where R± are 

= r ± (x>x )-T ± (x<x ) . (2.23) 
T±(x> x ) + T±(x <x ) v ; 

x is a dimensionless variable proportional to the energy, and xq is any fixed value in the 
energy interval. 

3. The if* -> W^b' process 

Following the formalism of Section 2 for the differential partial decay rate of the process 
(1.1) in the rest frame of H ± , we obtain 



dr± = r H± r° b 



1 ± a b rrif 



^L» vf . (3.1) 



(ptpv) 

T H ± is the partial decay width of the process H ± — > tb, assuming CP-violation in its vertex 

T H ±=C H (T inv ±T CP ), (3.2) 

where 

3a LO \ 1/2 (m 2 H ,m 2 t ,m 2 b ) g 2 

Ch = 4^P4 ' a - = ' (3 ' 3) 

r*"" = {Vt + vbiPtPb) ~ 2rn t m b y t y b , (3.4) 

T CP = [y t Re(5Y t CP ) + (y b Re(5Y b CP )] (p tPb ) - m t m b [y t Re(5Y b CP ) + y b Re(5Y t CP )] , (3.5) 

X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz, \ l ^ 2 {ra 2 H , m 2 ,m 2 ) m m 2 H — m 2 , (3.6) 

(PtP b ) = \( m2 H -m 2 - ml) ~ -(m 2 H - m 2 t ) , (3.7) 



,o = 9 2 {m 2 -m 2 w ){m 2 + 2m 2 w ) K ~ m 2 w )d cos 6 b ,; b , 

b 8m 2 w E t ' b '' b ' 16^1(1 -fa cos 6 , ; - 6 ,) 2 ' 



3.1 Angular distributions 

For the angular distributions of &'(&') when the 3-momentum of the t(t)-quark is along the 
z-axis one gets 



dr± = a 

dcos6 b , b , (1 — Pt COS 9y y) 2 



b _ C b J r i nv _| r C'P 



T mv ± + (3.9) 



£*(1 -PtCosOyy) m 2 t 
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where 

r ^l\Vt\{m\ - m 2 w ) 2 {m 2 t + 2m 2 w ) 

Mm 2 H m± w E 2 m t T t ' ^ 

a _ Wt\ I-., _ X 1/2 {m 2 H ,m 2 ,m 2 ) m 2 H - m 2 

~ ~E~t' H " 2^ -~2^-' (3 - 11} 

^ = + m 2 - m 2 ^ m 2 ^ + m 2 ^_ = + m 2 - m 2 ^ m 2 H - m 2 
1 2m H 2m H ' b 2m H 2m H ' 

We are interested in the CP-violating contributions to the loop corrections of the H^bt 
vertex (2.6). The quantities SY^ P and SY b CP enter the two independent combinations 
and V CP .One therefore needs two measurements to determine them. The decay rate 
asymmetry 5 b p for process (1.1) measures T CP given in (3.5), 

CP _ N V - Ny _ 

d b - Nb , +N - b , - V inv > V A6 > 

where N b ,t b ,\ are the total number of b' (V) in — ► bb'W^ decay. 

The CP-violating angular asymmetry measures the other combination V CP 

given in (2.17). We have 

The FB asymmetries are given by 

= A + ^™** \J b + J^ f^TfUP 1 ■ (3-15) 

Using the expansion 

V : —— I = V inv ± P CP + higher orders , (3.16) 

we get at one-loop level 

/ 2 2\ 

A b± - fi t + a h m t m H — - r — (V ±V ). (3.17) 

\ rn H m t) 

3.2 Energy distributions 

We write the energy distribution of b'(b') as a function of x = Ey /run (x = Eyjmji): 

dx 

where 



C E [c$ + cfx], (3.18) 



Ce ~ 2^ w m 2 H m t T t ' C ° " C ° ± C ° ' ^ ~ ° x ± ° l ' (3 ' 19) 

c inv = T inv _ ab ( m H+ m t) T invpinv^ (3 2Q ) 
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c[ nv = 2a b , ?^ m * 9 r inv P inv , cf p = 2a b -^i r -(T inv V CP + V CP V inv ). (3-22) 
{mi-m 2 w ) (mf-m^) 

The asymmetry AR CP also measures V ■ We choose xo = {x m in + x max )/2. Insert- 
ing the one-loop result 

R b± = \a h {m\ - m\) (V mv ± P CP ) (3.23) 

into eq.(2.23) gives 

AR° P = \a h {ml - m\)V CP . (3.24) 
4. The f/± ■"-v process 

In order to obtain the differential partial decay rate of (1.2), we fix the coordinate system 
such that the t-quark points in the direction of the z-axis and the 3-momenta of t and I 
determine the yz-plane: 

Pt,t = Wt,t\ (0,0,1), pi± = |^±|(0,sinfy±,cosfy±), 

Pb',b> = |Pfo',5'l( sin ^b'^' cos sin 6» 6 , 5 , sin 0^^/, cos 6» 6 , s , ) . (4.1) 
The angular distributions of are then given by 

(ft,tPi± 



dvf = v% r° 



1 ± aimt- 



(ptPi) J m t r t 



' '/'!>/ • (4.2) 



where 



rQ = ^K 2 - 2 (Pg')]^) , 6( P > w -m> w )d$ l± = JL ff , (4.3) 

1 2£ t m w IV lW w; ' (2tt) 4 4(m? - m^)m^ V ; 

{ Pt pi) = E t E l {l-p t <xxOi), (4.4) 

p = ^ „ nH p = rn 2 t -m 2 w 

1 2[E t {l-p\< X *0 l )-E v {l-< X *6 w )] b ' 2E t (l-p t coB0 v )' K ' } 

The angle Q\y is between p\ and jty: 

cos 9 ly = sin 0; sin 9y sin 0^/ + cos Q\ cos 6^/ . (4.6) 
The angular distributions of then reads 



dTf = (h_ 

dcos6i± (1 — PtcosO^) 2 



l _ W I r inv I -pCP 



^ „2 [r ffli)T)itM) I / r m»T,CP , -rCP<r,inv\-\ ( E b(^ + cos fy±) (g*gb) A 

a/m t 1 P ± (1 /■> +1 A> J — — — — 

\E t (l - p t cos9i±) mf J 
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where 



C t = -alm w \p t \(m 2 t - m 2 w f x (4.8) 
[6(1 - PlfEj + 3mj - hmWw + ~ 3(1 - 0?)E?(3m? t - 2m 2 w )} 
^m\m t Y t Y w El\m\ - m 2 w - (1 - P?)E?}* 

Tw is the total decay width of the W boson. 

As there is no CP violation in t — ► blv decay, the decay rate asymmetry 5p p for process 
(1.2) will measure the same combination T CP : 

' N l+ + N t - T mv v ; 

where iVj± are the total number of in i7 ± — ► bb'l^v decay. 

For the CP-violating FB asymmetry AA CP of the process (1.2) we obtain 

AAf p = 2a im 2 t rnl ( m }~ m l v CP , (4.10) 
' (m^ + m 2 ) 2 

and the FB asymmetries are at one-loop level 

(2 2\ 

AFB _ a , ^ ^2 2 \ m H ~ m t j taint , -pCP^ /* ii\ 

A l± - fa + aim t m H + 2 (P ±P ). (4.11) 

Notice, that the only difference between the expressions (3.14) and (4.10) is the coefficient 
a b in (3.14) and a\ in (4.10). These coefficients are only connected to the polarization of 
the t quark. Because of the fact that a b = 0.38 and a\ = 1, one would expect a bigger 
effect measuring (4.10). 

5. Numerical Results 

Here we present a numerical analysis of the discussed asymmetries. First we analyze the 
CP-violating asymmetries 5^f , AA^f and AR^ P . Further, we study the FB asymmetries 
A£f ± and R b ± needed for AAgf and AR^ P . 

5.1 The CP- violating asymmetries 

The expressions V CP andr CP , (2.17) and (3.5), are linear combinations of the CP violating 
form factors Ke(5Y t CP ) and Re(<5Y" 6 ). Therefore, we need to measure: 1) the decay 
rate asymmetries S^f which are proportional to T CP and 2) the angular and/or energy 
asymmetries which are proportional to V CP . 

As there is no CP violation in t — ► bW, the decay rate asymmetries for (1.1) and (1.2) 
are equal to the decay rate asymmetry for — > tb. We denote it, following ref.[5], by 
5 CP : 

5 CP = 5<; p = 5f p = ?Z. (5.1) 
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The angular and energy asymmetries are not independent either. As seen from (3.14) 
and (4.10), the angular asymmetries for leptons and 6-quarks are related by: 

AA cr = AA£ P « 2.6 AA% P . (5.2) 

Further, (3.14) and (3.24) lead to a simple relation between the 6-quark energy and angular 
asymmetries: 

A^^ ^Vi^ ^r, (5.3) 
4 m H+ 

which implies that for m H + > mt, AR^ P is bigger than AA^ P . Thus, in general, AR^ P 
is the biggest asymmetry of 2) for m H + > 490 GeV. Fig. 2 illustrates the relative size of 
the asymmetries AA?F and AR^ P as a function of m H + Note that the relations (5.2) and 
(5.3) are independent of tan f). 




400 600 800 1000 1200 1400 

m H + [GeV] 

Figure 2: The ratio AR^ P /AA^ P as a function of m H + 

Therefore, in the following we shall discuss only the decay rate asymmetry 8 and the 
energy asymmetry AR^ P . (Because of a conjugation error in our paper [5] in the formula 
for the tbg vertex, see the Appendix B, we have redone the numerical analysis for 8 CP .) 
The purpose of our analysis is to determine the size of the asymmetries as functions of 
rrifj+ and tan (3, being the most important parameters of the Higgs sector in MSSM, for 
different values of the CP-violating phases. 

The sources of CP violation in our processes are the one-loop corrections to the H + tb 
vertex with intermediate SUSY particles, see Fig. la, lb, Id, le of [5] and the self-energy 
graph with fv T . (The corrections due to Fig. lc and Fig. If of [5] are of higher order and we 
do not consider them here.) In order not to deal with too many phases we assume the GUT 
relation between M\ and M 2 which fixes <$>M\ = 0. According to the experimental limits 
on the electric and neutron EDM's, we take ^ = or ^ = 7r/10. Thus, the remaining 
CP-violating phases in our study are the phases of A t , A^ and A T which we shall vary. If 
not specified otherwise, we fix the following values for the other MSSM parameters: 

M 2 = 300 GeV, M 3 = 745 GeV, M (j = M 6 = M 3 = M E = M L = 350 GeV, 
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172 


519 
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464 
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402 


344 



Table 1: Masses of the sparticles (in GeV) for the parameter set (5.4) together with cf>A t = 4>A b = 
tt/2 and M = 0. 



H = -700 GeV, \A t \ = \A h \ = \A T \ = 700 GeV (5.4) 

The relevant masses of the sparticles for the choice (5.4) and tan (3 = 5 or 30 are given 
in Table 1. For the case with = 7r/10 and the other parameters unchanged, all masses 
do not change by more than 1 GeV from those given in Table 1, except for mi = 187 GeV 
and mr = 515 GeV for tan/3 = 5, and = 176 GeV for tan = 30. Note that ^ = 7r/10 
implies fi = — 700e J7r / 10 GeV. We have used running top and bottom Yukawa couplings, 
calculated at the scale Q = m H + , see Appendix A. We have checked that the asymmetries 
have only a very weak dependence on the scale Q. 

Fig. 3 shows the asymmetries 5 and Ai?^ p as functions of m H + for <pA t = tt/2, 
4>A b = and (f>n = 0. As one can see, for tan/3 = 5 the decay rate asymmetry 8 CP goes up 
to 20%, while ARg p reaches 8% for the same values of the parameters. The asymmetries 
strongly depend on tan/3 and they quickly decrease as tan/3 increases. Our studies have 




Ai?r 0.08 




400 600 800 



1000 1200 1400 

m H + [GeV] 



400 600 800 1000 1200 1400 

m H + [GeV] 



Figure 3: The asymmetries S CP and AR^ P as a function of m H + for <f>A t 
The red, blue, and green lines are for tan/3 = 5, 10, and 30. 



tt/2, ^ 



= ^ = 0. 



shown that the most important CP-violating phase is (j)^ 
pendence on <pA b and <j)A T - We therefore take them zero. 



There is only a very weak de- 



The main contributions to both 8 and AR^ P come from the self-energy graph with 
stop-sbottom. The vertex graph with stop-sbottom-gluino also gives a non-zero contribu- 
tion. Their contributions are shown in Fig. 4. The contribution of the rest of the graphs is 
negligible. This justifies the use of the GUT relation which fixes <f>Mi = 0, and it explains 
the weak dependence on <f> Ar . 
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b d 

400 600 800 1000 1200 1400 400 600 800 1000 1200 1400 



m H + [GeV] m H + [GeV] 

Figure 4: The contribution of the tb self-energy (red line) , tbg vertex contribution (blue line) and 
the sum of the other (green line) diagrams to S CP and Ai?£" p as a function of m H + for tan j3 = 5 
and 4> At = 7t/2, (j) Ab = fa = 0. 



Up to now, all the analyses are done for M3 = m<j = 744 GeV. In Fig. 5, we show the 
dependence of 5 on the gluino mass. In general, 5 CP gets small with increasing gluino 
mass. 




0.14 u — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — u 
300 500 700 900 

nig [GeV] 



Figure 5: 5 CP as a function of the gluino mass for tan/3 = 5. The red line is for m H + = 600 GeV, 
the blue line is for m H + = 800 GeV and the green line is for m#+ = 1000 GeV. 

Let us now allow a non-zero phase of [i. We take a very small phase, <^>« = 7r/10, in 
order not to be in contradiction with the experimental data. As can be seen in Fig. 6, the 
asymmetries can increase up to 25% for 5 CP and 10% for AR% P . The discussed asymme- 
tries 5 CP and Aii^ p show a very strong dependence on the sign of \i. As noted above, 
our analysis is done for /i = —700 (see (5.4)), however if u, changes sign, \x = 700, all 
asymmetries become extremely small. 

5.2 The P-violating asymmetries 

When discussing the possibilities to measure AA^f and AR^ P , it is also important to 
know the size of the FB asymmetries A F P,, (3.17) and (4.11), and of the energy asymme- 
try Rb±, (3.23), that enter the corresponding CP-violating asymmetries. 
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-U.Uki 

400 600 800 1000 1200 1400 4 0° 600 800 1000 1200 1400 



m H + [GeV] m H + [GeV] 

Figure 6: The asymmetries S CP and Ai?p p as a function of m H + for <pA t = tt/2, </u b = and a 
non zero phase of /x, </> M = tt/10. The red, blue and green lines are for tan/3 = 5, 10, and 30. 



Afrp± and Rb± are determined by the polarization of the t-quark in — * tb 
decays. As the Lagrangian violates parity, these asymmetries appear already at tree level 
and thus should be rather large. 

Neglecting the loop induced CP-violating part T> in (3.17), (4.11), and (3.23), we 

get 

Ainv_i( A FB, A FB\_a, 2 2 \ m H+ ~ m t ) v inv (r. r\ 

A b,i - r + A b,i-) - Pt + aum t m H+ —2 y , (D.dj 

z \ ' n H+ ~ l ~ t ) 

i?r = \ (R b+ + = \a b (m 2 H+ ~ m 2 t )V™ ~ ^4^4 ■ ( 5 - 6 ) 

2 4 ^ Vt+Vb 

Thus AAgf , (3.14, 4.10), and Ai?f p , (4.10), are determined by V CP , while Ajf and 
are determined by J> mv ^ and there is no a priori reason to expect that their tan/3 and 
m H + dependences will be the same. A % ™ and R 1 ™ are tree-level quantities. Including the 
one-loop corrections to these would require the full renormalization of the process which is 
beyond the scope of this paper. 

In Fig. 7 we present A™ v and A 1 ™ as a function of m H + for tan/3 = 5,10, and 30, and 
in Fig. 8 we show R l b nv as a function of tan/3. Note that at tree-level AT? depends only 
on tan /3 and m H + , and only on tan /3. 

6. Conclusions 

We have calculated the CP-violating decay rate, forward-backward and energy asymmetries 
between H + -> It -> bb' W+(-^ bb'l+Ui) and JET" -> bt -> bb'W'^ bl'l'vi). They are 
induced by loop corrections in the H^tb- vertex. The CP violating forward-backward and 
energy asymmetries are determined by the polarization of the top quark and are therefore 
related. We have shown that it is necessary to measure both the decay rate asymmetry 
5 CP and the forward-backward or the energy asymmetry to get the maximal information 
on the CP-violating parts of the decay amplitude. We have performed a detailed numerical 
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Figure 7: The forward-backward asymmetries A™ v and A % ™ as a function of m H + for tan/3 = 
5 (red), 10 (blue), and 30 (green). 
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Figure 8: The energy asymmetry R l b nv as a function of tan j3 

analysis of these quantities. An important improvement with running Yukawa couplings 
at the m H + scale has been made. The asymmetries are most sensitive to the phase <pA t - 
The asymmetries reach their maximum for tan/3 = 5 and \x = —700 GeV. The decay rate 
asymmetry can go up to 25%, the others up to 10%. The main contribution comes from 
the self-energy diagram with stop and sbottom exchange. We have also calculated the 
P-violating asymmetries at tree level. 

We want to add a few remarks on the measurability of these asymmetries. In principle, 
the production rate for H at LHC is not so small being 0.2 pb for = 500 GeV and 
tan/3 = 30 [10, 11]. The main production process is due to g b — > H + i. Because of the 
large background, the actual signal production rate is strongly reduced. According to [10], 
one can expect N = 733 signals with N/\fB = 12.6 for = 500 GeV, tan /3 = 50 for a 
luminosity C = 100 fb _1 . The statistical significance y/~N A to measure an asymmetry A 
of several percent might be too low for a clear observation of CP violation in H + decays 
at LHC in the first stage. However, at SLHC for which a luminosity of 1000-3000 fb _1 
is designed, such a measurement would be worth of being performed. Here we have only 
considered CP violation in the H + decays. However, similar graphs are also present in 
the production process g b —* H + t [8]. One would expect a CP-violating asymmetry 
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of similar size. The total asymmetry in production and decay would be approximately 
additive, A tatm = A prod , + ^decay- 
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A. Running Yukawa couplings 



For clarity, we present all formulas used for programming the running top and bottom 
Yukawa couplings, h b and ht, respectively. The Lagrangian for the H ± tb interactions reads 



C Hqq = H + t{y* b P R + y t P L ) b + H-b (y* t P R + y b P L ) * + ..., 
with the DR running top and bottom Yukawa couplings in the MSSM, 

y b = h b sin , y t = h t cos (3 , 
given at the scale Q = m H + in our studied case, and 



h b 



\/2mw cos P ' 



\^2mw sin (3 



(A.l) 



(A.2) 



(A.3) 



In [4] it is shown that within an effective theory approach large finite scale independent 
parts can be resummed, which in case of complex MSSM input parameters leads to complex 
h b and h t . Effective means that the masses of the particles in the loops are much bigger 
than those of the in- and outgoing particles so that these states can be integrated out in the 
Lagrangian. In our case, we are interested in additional open channels, e. g. H + — > tb. This 
implies that the resummation is not applicable here. But we can improve our calculation 
by using full one-loop running quark masses with some higher order improvements of the 
gluonic part. Note, that m^ R (Q) can always be made real by field redefinition [4] and 
therefore also h q is real in our case. 

We take as input set the bottom mass m^ s (mfe) = 4.2 GeV, the mass for the top 
quark is the pole mass, m]? ole = 171.4 GeV, the strong coupling is a^ s (mz) = 0.1176, 
m z = 91.1876 GeV, and m w = 80.406 GeV 

First we want to have the DR bottom mass of the Standard Model at the scale Q (see 
for comparison eq. (26) in [12]), 



<f M (Q) = m^(Q) 



MS, 



a 



DR 



3vr 



23a,' DR 
72vr 



(A.4) 



with mf ls (Q) = m b (Q)sM given in [13, 14]. Adding the loop contributions due to supersym- 
metric and heavy SM particles, denoted by Am^ xtra (calculated in the DR renormalization 
scheme), we get the full one- loop DR running bottom mass (with some higher-order im- 
provements) within the MSSM, 



m DK < 
m b 



(Q) = mgf M (Q) + Amr(Q) 



The DR running top mass we get from 



(Q) 



m 



pole 



Am, 

1 + 



(i) 



Am 



(2,3) 



m t 



(A.5) 



(A.6) 
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where Am^ is the full one-loop contribution to mt (calculated in the DR renormalization 
scheme) and Amf' 9 ' is the gluonic two-loop contribution, 



Am 



(2,9) 



a s (Q)y {8ir 2 2011 16tt 2 log(2) 8£(3) 



4-7T 



18 



+ 



82L + 22L 2 , (A.7) 



with L = log(Q 2 /m 2 ), see [15]. 



B. Squark-quark gluino contribution 

In Appendix B of [5], eq. (62) is incorrect and therefore also eqs. (14, 15). For clarification, 
the definition of the squark rotation matrix is essential. In this work and in [5], one 
has 

q a = R*i a qi with a = L,R and i = 1,2. (B.l) 
Hence, the squark-quark-gluino interaction ( eq. (62) of [5] ) is given by 



l q¥g = -V2g s T«[?(Rt*e-^P L -Rt*e L ^P R ) qs ql 

+ q s (Rl ~ 4i e-^ 3 P L ) ~g a qi,t 



(B.2) 



The contribution from the diagram with a stop, a sbottom, and a gluino (in [5] 
eqs. (14,15)) is 



4 a, 



Re5Y b CP {tibfg) = --^ [m- g Im^i^i?!* e^) Im(C ) 

O 7T 

+m t ImiGajR^Rl*) Im(Ci) + m b Im^i?^*) Im(C 2 ) ] , (B.3) 



4 a, 



ReSY^iUb.g) = --— [m g Im(G Atj RX* e~ 1 ^) Im(C ) 

+m t ImiG^jRliRf*) Im(Ci) + m b Im^i^l,*) Im(C 2 ) ] , (B.4) 



with C x = C x (m 2 , m H+ , m l > m § > m | » m f . ) • 
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